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j> ■ Abstract 

O 

We investigate the current status of the long-wavelength vacuum oscillation 
i^) solution to the solar neutrino problem and to what extent the presence of 

CD ■ a third neutrino can affect and modify it. Assuming that the smaller mass 

squared difference that can induce such oscillations, Am 2 2 , is in the range 
10 -11 — 1CP 8 eV 2 and the larger one, Am^, in the range relevant to atmo- 
spheric neutrino observations, we analyze the most recent solar neutrino data 
coming from Homestake, SAGE, GALLEX, GNO and Super-Kamiokande ex- 
*^ periments in the context of three neutrino generations. We include in our 

vacuum oscillation analysis the MSW effect in the Sun, which is relevant for 
^ . some of the parameter space scrutinized. We have also performed, as an 

extreme exercise, the fit without Homestake data. While we found that the 
MSW effect basically does not affect the best fitted parameters, it significantly 
modifies the allowed parameter space for Am 2 2 larger than ~ 3 x 10 -10 eV 2 , 
in good agreement with the result obtained by A. Friedland in the case of two 
generations. Although the presence of a third neutrino does not essentially 
improve the quality of the fit, the solar neutrino data alone can give an upper 
bound on #13, which is constrained to be less than ~ 60° at 95 % C.L. 
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I. INTRODUCTION 



The solar neutrino problem (SNP) [|IJ seems now to be established as a definite sig- 
nal of non-standard neutrino properties by the four first generation solar neutrino experi- 
ments, Homestake 0, Kamiokande [||], GALLEX [|], SAGE and by the new generation 
higher statistics solar neutrino experiments, Super-Kamiokande || and GNO @, which have 
strengthened the existence of the SNP. Astrophysical explanations of the SNP, which requires 
significant deviation from the standard solar model (SSM) ||||, are highly excluded with 
the current solar neutrino data ||10||. Moreover, there is an excellent agreement between the 



sound velocity predicted by the SSM and that obtained from the recent helioseismological 
observations |TT|, which supports the rigidness of the SSM. 

Albeit electron neutrinos, produced in the Sun, most certainly vanish on their way to 
the Earth, which is the dynamical origin of the process that promotes their disappearance 
is yet to be completely clarified. It has been discussed that the SNP can be nicely explained 
by the simplest extension of the standard electroweak model which invokes neutrino mass 
and flavor mixing fll2" |. The most plausible solutions in this context are provided either by 
the matter enhanced resonant neutrino conversion, the MSW effect [ 13| , or by the vacuum 
oscillation [II] with a typical wavelength as long as the Sun-Earth distance [|15 . Recent 
detailed analyses and discussions of the solar neutrino data based on these mechanisms, i.e. 
the MSW and long-wavelength vacuum oscillation (LVO) solutions, in the context of two 
neutrino generations, can be found, in Refs. [|16|,|l7| and Refs. fL6Ul8f , respectively. Earlier 
detailed analyses, prior to the Super-Kamiokande experimental results, can be found in 
Refs. [ 1^,20]. Discussions on other possibilities to explain the SNP by invoking more exotic 
properties such as neutrino magnetic moment, flavor changing interactions or violation of 



the equivalence principle, can be found, for example, in Ref. |21 

Although it seems at present that pure z/ M — > v T oscillations in vacuum are quite enough to 
account for the atmospheric neutrino anomaly |22| . |23| , the possibility of having contributions 
from non-negligible v^ — » v e oscillations is still not discarded |24 ]25|j , even after taking into 
account the constraints coming from the CHOOZ reactor experiment |p5 ]. Moreover, the 
existence of at least three neutrino flavors, which is one of the most impressive results from 
the LEP experiments |2"7|,2"5], makes it unavoidable to try to understand neutrino oscillations 
in a full three generation scenario. 

In the framework of three generations of massive neutrinos subjected to flavor mixing, 
one has to deal, in general, with six variables to study the neutrino oscillation phenomena, 
namely, three mixing angles, #12, #23, #13 and one CP violating phase 5, which relate mass 
and flavor eigenstates, and two (independent) mass squared differences, which can be chosen 
as Am 2 2 = m\ — m\ and Am^ = m\ — m\. If we assume that the smaller mass squared 
difference, defined as Amf 2 , is in the range relevant to solar neutrino oscillation, i.e. Am 2 2 ~ 
1CT 11 — 1CT 8 eV 2 in the case of the LVO solution or Am 2 2 ~ 1CT 8 — 1CT 5 eV 2 in the case 
of the MSW solution, and the larger one, Am 23 , is in the range relevant to the atmospheric 



neutrino observations, 



i.e. 



Am 2 3 



10" 



10 eV , only three of these variables become 



significant in practice to the solar neutrino investigation: #12, #13 and Am 2 2 [f29"H 3~T 



Under such assumption, a detailed analysis of the three flavor LVO solution to the SNP 
was performed in Ref. |3^] in which the best fit was obtained when # 13 is zero, this angle 
being constrained to be less than ~ 30° at 95 % C.L. by the solar neutrino data alone. The 



authors of Ref. |32| have in addition taken into account the atmospheric neutrino results, and 
found that the combined three generation fit does not lead to an allowed region appreciably 
different from the one obtained by performing two separate effective two-neutrino fits to the 



solar and atmospheric neutrino data. See, for instance, Ref. |H|] for earlier discussions on 
the three flavor LVO solution to the SNP. 

On the other hand, a detailed three flavor analysis of the MSW solution to the SNP was 
performed in Ref. JS3| and not long ago its updated version has been reported in Ref. [|35| . 
The best fit obtained by this study permitted to constrain # 13 < 55° — 60° at 95 % C.L. by 
the solar neutrino data on its own, without taking into account the atmospheric neutrino 
observations p2| or the CHOOZ result [p6]. More recently, a complete analysis considering 



four neutrinos, three active and one sterile, was performed in the context of both MSW as 



well as LVO oscillation solutions [Q, and bounds on mixing angles were again obtained 
solely by the solar neutrino data. 

In this same spirit we re-examine and update in this paper the long-wavelength vacuum 
oscillation solution to the SNP in the two generation as well as in the three generation 
frameworks. In this study, we do not take into account the constraints on #13 coming from 
the Super-Kamiokande atmospheric neutrino observations [^2| nor by the CHOOZ reactor 



experiment |26|| , whose combination gives an upper bound on #13 corresponding to sin 6 



1.3 



< 



few % [p4p5| . The main idea of this paper is to constrain # 13 by LVO and the solar neutrino 
data alone, as it has already been done for the atmospheric neutrino pi| , |2"5| as well as for 
the MSW solar neutrino solutions [35 . 

Moreover, we investigate the effect of taking the 8 B neutrino flux as known only up to 
a normalization factor, /b, renormalizing, in this case, the 7 Be flux accordingly, assuming 
the solar temperature power-law [37|. Finally, we also perform an extreme check analysis 



by neglecting completely the solar neutrino rate measured by the Homestake experiment, 
as considered, for instance, in Ref. 



35 



on the account of it being the only radiochemical 
experiment which has not been calibrated with a radioactive source. 

There are new ingredients we include in this work, which were not considered in previous 
analyses of the three flavor LVO solution. We extend the analysis to mass squared differences 
up to 10~ 8 eV 2 in order to cover the whole relevant LVO parameter space. Above this value 
the oscillation probability of solar neutrinos essentially does not depend on the distance 
between the Sun and the Earth as discussed in Ref. [j39fl , and we do not regard such case 
as a LVO solution to the SNP. Thereby, we include in our estimations the MSW effect in 
the Sun in the context of the LVO solution, whose importance was first pointed out in 
Ref. g0| for Aml 2 /E larger than 3 x lCT 10 eV 2 /MeV. The relevance of the MSW effect 



in determining the LVO allowed parameter space has been recently discussed in Ref ET 
Due to the presence of matter effect, we cover the whole relevant range of mixing angle 
< 612 < 7r/2 ['11] in contrast to the range usually considered (0 < 612 < 7r/4) for the LVO 
solution. 

We have used in our analysis the most recent solar neutrino data from the five solar 
neutrino experiments^ : the total rates measured by the Homestake chlorine (CI) detec- 
tor @, by the GALLEX @], GNO and SAGE § gallium (Ga) detectors, which are 
summarized in Table |, as well as all the observations performed by the high-statistics Super- 
Kamiokande (SK) water Cherenkov detector (total rate, energy spectrum and zenith angle 
dependence) H. The pioneering Homestake experiment which has its energy threshold at 



1 For the sake of simplicity, we have not included the data from the Kamiokande experiment ||, 
which can be justified by the fact that its result is consistent with that of SK and its errors are 
much larger than that of the SK experiment. 



Ey = 0.814 MeV |l] is mostly sensitive to 8 B neutrinos (~ 76 % of the total contributions), 
and also can detect 7 Be neutrinos (~ 15 %). The gallium detectors have their threshold 
at E v = 0.233 MeV [| and are sensitive to pp (~ 54 %), 7 Be (~ 27 %), and also to 8 B 
neutrinos (~ 10 %). Super-Kamiokande which has its energy threshold at E e = 5.5 MeV, 
where E e is the total energy of the recoil electron, is only sensitive to 8 B neutrinos^]. 

This paper is organized as follows. In Sec. [TT] we present the oscillation formalism we will 
use to pave the way to the two and three neutrino generation analyses of the data. In Sec. 
I IT] we present and discuss the general features of our data analysis. In Sec. [TV] we revise and 
update the two-generation LVO solution also including the MSW effect in the Sun. In Sec. 
|V] we extend the analysis to three-generation neutrino oscillation taking into account the 
MSW effect. In Sec. |V]]we examine the seasonal effect that can be expected at SK, Borexino 
and KamLAND experiments for the best fitted parameters of the LVO solutions we found. 
In Sec. |V11| we discuss the LVO solution with two and three neutrino flavors disregarding 



the chlorine data. Finally we draw our conclusions in Sec. [VIII . 



II. OSCILLATION FORMALISM 



In this section we describe the framework we will use in the scope of this paper. 



A. Two Flavor Case 

For the two generation case in which v e is mixed only with v^ (the same argument holds 
also for the case of v e — v T mixing), the evolution equation of neutrino traveling through 
matter can be written as 



. d 

dx 



V e (x) ±A 12 sin2fl 12 
|A 12 sin 2#i2 Ai 2 cos 2#i2 



(1) 



where Ai 2 = Aml 2 /2E, with the mass squared difference of the two neutrino mass eigen- 
states Atoj 2 = m 2 — m\ , E is the neutrino energy, #i 2 is the mixing angle, 



V e {x) = V2G F N e (x) ~ 7.6 x 10" 14 



N P 



1 mol/cc 



eV 2 



(2) 



is the matter potential for v e with Gf and N e being the Fermi constant and the electron 
number density, respectively. 

This equation, in the case of vacuum (V e = 0), simplifies to lead to the well known 
formula for the v e survival probability, 



2g,vac^e 



sin 2 26*12 sin 2 



•A 



12 



L 



1 - sin 2 2#i 2 sin 2 1.27 
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[ eV 2 _ 




"MeV 

[ E 




L- 
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(3) 



2 In this work, we do not take into account the possible enhancement of hep neutrinos [42] in the 
SK spectrum analysis since the current SK data do not seem to conclude on the necessity of such 
enhancement [pi. 



where L is the distance traveled by the neutrino. It is clear that the probability P2gvac is 



invariant under the transformation Am 2 2 ■ 
in vacuum it is sufficient to assume Am 2 2 



-Amf 2 as well as #12 



:g,vac 

| — 9 12 and therefore, 



> and < ^12 < f to account for all physical 
situations. 

From Eq. (|3|) we can estimate the vacuum oscillation length as, 



4irE 
Am 2 12 ~~ 

0.165 L 



2.47 x 10 10 



E 
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Ll MeV 



10 



1 MeV 
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12 



10- lu eV 2 
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(4) 



where Lo = 1 AU ~ 1.496 x 10 m, is the astronomical unit, namely, the mean Earth-Sun 
distance. 

If matter is present, the physics described by the neutrino evolution given in Eq. ([!]) 



is not invariant under either transformation AmJ 2 



Am? 9 or 9 



'12 



12 



TV Q 

2 ~ y l2- 



However, 

we can show that the physical consequence of Eq. (|l|) is invariant under the transformation 
6?i2 — ^ 7r — #12 as well as (Am 2 2 ,6>i 2 ) — ► (-Am^, | — #i 2 ). Therefore, in order to cover all 



the possible physically meaningful parameter region we can assume [31.35 



Am| 2 > 0, 



IX 



< 0i a < - 



(5) 



This is the range we consider in this work. 

In principle, neutrinos on their way from the interior of the Sun to the detector on the 
Earth can be influenced by the solar as well as by the Earth matter potentials. 

Let us initially consider what can happen in the Sun. In Fig. |l|(a) we show the iso- 
survival probability contours of v e at the solar surface as a function of Am 2 2 / ' E and sin 2 #i 2 
for pure vacuum oscillation, ignoring the matter effect, assuming that neutrinos are created 
in the solar center. In Fig. |l](b) we show the same type of contours but taking into account 
the matter effect. The probabilities were obtained by numerically integrating Eq. (JI|), using 
the electron number density one can find in Ref. f43|, again assuming that neutrinos are 
created in the center of the Sun. 

From these two plots, we immediately see that the probability at the solar surface can 
be significantly different for these two cases if Am\ 2 /E > few xlO -10 eV 2 /MeV. Therefore, 
one can not disregard the matter influence in the v e survival probability calculation for 
values of Am 2 2 / 'E in this range. For smaller values of Am\ 2 /E the oscillation effect is small 
and the survival probabilities are ~ 1, in both cases. For Am\ 2 /E ~ (9(10 -9 ) eV 2 /MeV 
or smaller, the role of the solar matter is to suppress oscillations, as Am\ 2 / E increases 
the MSW resonance effect comes into play and strong conversion can occur. We note that 
for values of Am 2 2 > 10~ 9 eV 2 , the MSW conversion can be significant for lower energy 
neutrinos, such as pp and 7 Be ones, but it will be very weak for most of the 8 B neutrinos. In 
fact, we can see from Fig. [I] (b), that even for Am 2 2 = 10~ 8 eV 2 , the highest value considered 
in this work, if we take E ~ 10 MeV, the typical 8 B neutrino energy relevant for SK, the 
survival probability of u e at the solar surface is very close to one. This implies that at such 
high energies, neutrinos exit from the Sun essentially as v e . 

Let us now examine what happens in the Earth. For the averaged electron number 
density found in the Earth's mantle (N e ~ 2 mol/cc) or core (iV e ~ 5 mol/cc), the matter 
potential V e is always much larger than A12, unless Am 2 2 is very close to 10 -8 eV 2 for 
the lowest energy (i.e. pp) neutrinos (see Eq. (§)), leading to a strong suppression of the 



effective mixing angle in the Earth matter. This implies that no appreciable v. 



/(,r 



v e or 



v e — ¥ u ^,t) re-generation effect can occur in the Earth, since this can promote, at the most, a 
few % change in the probability for a very limited range of mixing parameters (A12 ~ 10~ 7 
eV 2 /MeV and large 612) ■ Hence, in this work we only take into account the MSW effect in 
the Sun and neglect it completely in the Earth. 

After accounting for the MSW effect in the Sun, the v e survival probability at the Earth 
given in Eq. 



is modified as follows 40 



P. 



2g 



V„ 



u e )= cos 2 9 12 \A U1 (R Q )\ 2 + sm 2 9 12 \A I/2 (R Q )\ 2 

+ sm29 12 \A Ul (R Q )A l/2 (R Q )\cos(A 12 L + (3) 



(6) 



where A Ul (R & ) and A U2 (R & ) are respectively the amplitudes of the neutrino state to be found 
in the mass eigenstate z/i and v 2 at the solar surface, L is the distance between the solar 
surface and the detection point, and (3 = Aig[A Ul (R & )Al 2 (R Q )] corresponds to the phase 
developed between v\ and v 2 after the neutrinos pass the resonance point until they reach 
the solar surface. We remark that for the values of AmJ 2 we are interested in this work, in 
the region where z/ e 's are created, typically r < 0.3R Q , because of the condition V e 3> Ai 2 , 
they essentially coincide with the mass eigenstate z/i(~ u e ), and no oscillation occurs within 
this region. Hence the final probability amplitudes A Ul (R Q ) and A U2 (R Q ) as well as (3 do not 
depend on the exact production point, allowing us to assume that all neutrinos are created 
in the solar center [ffOf . 

The probability in Eq. (pf) can be rewritten in terms of the flavor amplitudes as, 



P 2 gOe -» V e 



= \A Ve (R Q ] 
-sin29 12 \A Ue (R Q )A Ux (R G 



1- sin 2 2#i 9 sin 2 f -^L 



2,,;,.-2 0i ) ,.;,,:< I ^l 2 .^ 



+ |v4 1/x ( J R )| 2 sin"2^ 12 sin 



2 008 26*12 sin 



A 



12 



L ) cos f3' - sin ( A 12 L) sin f3' 



(7) 



where A Ue (R Q ) and A Ux (R Q ) are respectively the probability amplitudes of the neutrino to be 
found in the state v e and z/ x (x = /i, r) at the solar surface and f3' = Arg[y4y e (i? )y4*^(i? )], 
is the phase difference developed in the flavor basis corresponding to the one we have in Eq. 
(|) between the mass eigenstates, namely, (3. We note that the following relation holds, 



\A„ e (R 



©y 



cos 2 9 12 \A U1 (R Q )\ 2 + sin 2 9 12 \A U2 (R Q )\ 2 + sm29 12 \A Ul (R Q )A U2 (R Q )\cos[3. (8) 



In this work, we first compute the values of A Ue (R Q ) and A Ux (Rq) by numerically inte- 
grating the neutrino evolution equation in Eq. (|l|) with the solar electron density taken from 
Ref. ^3] and then compute the final probability at the Earth by using the expression in Eq. 
(^). We take into account the effect due to the eccentricity of the Earth orbit around the 
Sun by taking the time average over one year of the probability using the time dependent 
distance L(t), 



L(t) 



1 — e cos 2tt 



T 



+ 0(e 2 ), 



(9) 



where Lq = 1 AU, e = 0.0167 is the orbit eccentricity and T = one year. We neglect the 
correction terms of 0(e 2 ). 

Following Ref. [J4l| , we restrict the range of Am 2 from 10~ n to 10 -8 eV 2 . One can easily 
estimate from Eqs. ([|) and (f|), that for the typical energy of 8 B neutrinos, ~ 10 MeV, 



in order to have appreciable oscillation effects for solar neutrinos, the Am 2 2 value must be 



10 n eV 2 . This sets our lower limit on Am 2 2 . 



larger than of the order 

It is known (see, for e.g. , Ref. |41fl ) that if Am 2 2 is larger than ~ 10~ 8 eV 2 , the final v e 
survival probability at the Earth, averaged over the neutrino energy, would not practically 
depend on the precise value of the Sun-Earth distance since the variation in the probability 
due to the energy spread is large enough, even for 7 Be neutrinos, to average out the vacuum 
oscillation effect. Since we do not regard such a case as a LVO solution, this sets our upper 
limit on Am 2 2 . Also, to cover values of Am 2 2 up to 10 -8 eV 2 means to examine the effect 
of the solar matter on the LVO solution to the SNP whose relevance was already discussed 
in Refs. HIT]. 



B. Three Flavor Case 

For the case when we consider three flavor mixing, the evolution equation of neutrinos 
in matter can be written as, 
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(10) 



Am 



where Ay- e: ~ 2 ^ j and we take the representation of the Maki-Nakagawa-Sakata |L2| mixing 
matrix U which is the leptonic analog of the Cabbibo-Kobayashi-Maskawa matrix in the 
quark sector [04| as, 



TJ _ „i^7023 pi^50l3 „iM9l2 



(11) 



Explicitly 0, 



U 



C12C13 
-S12C23 - C12S23S13 

■S12-523 ~~ C12C23S13 



S12C13 S13 

C12C23 - S12S23S13 S 23Cl3 
— C12S23 — S12C23S13 C23C13 



where Aj are the SU(3) Gell-Mann's matrices and c. 



Am 2 3 
observations 123 



COS v%j , S{j 



(12) 



sindij. We assume 



to be in the range ~ 10 -3 — 10 -2 eV 2 , which is consistent with atmospheric neutrino 

Because of this assumption, the relation V e , A 12 -C A 2 3 ~ A 13 holds. 

In this case, A 13 is the dominant term in the Hamiltonian matrix and this leads to 

the decoupling of v% from the remaining two states. This means that the oscillation which 

is driven by the A 13 term can be simply averaged out in the final survival probability of 

electron neutrinos at the Earth, PzAy e -^ 



yielding the relation [p9| ,J3 



% 



v„ 



v„ 



sin 4 9 1 3 + cos 4 0i3 ■ P 2si 



Ve), 



(13) 



where i^g^e — > v e ) is defined as the survival probability of the following effective two 
generation system p9| , |3llj30| , 



3 Here we neglect the possible CP violating phase since this phase will not affect the v e survival 
probability, even if we take into account the next-to-leading order corrections in electroweak inter- 
actions 
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dx 
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where za, 



o: 



e, //, r) is defined as, 



v r 



e -JA 5 6»i3 e -'JA 7 6»23 



z/ e 
z^„ 



(15) 



We compute P 2g 



/A. 



P e ) by using the expression in Eq. (|7|) just by replacing the solar 
electron density as N e (r) — > cos 2 9±3 ■ N e {r). 

We note that we have three independent parameters Am 2 2 , 9\2 and #13 which can be 
fitted or constrained by the experimental data. Regarding #13, it is clear that the probability 
in Eq. ( p~3| ) is invariant under 13 — > it — 13 allowing us to restrict the range of 8 13 to 
< 613 < 7r/2. As for Am 2 2 and $12, we again consider the range given in Eq. ([50 for the 
case of two generations, since the same argument discussed in the previous subsection also 

holds for P2g(^e — ► v^. 



III. DATA ANALYSIS 

The general idea is to perform a x 2 analysis to fit the oscillation free parameters (812 and 
Am 2 2 , in the case of two generations, and 812, 813 and Am 2 2 in the case of three generations) 
with observed experimental data. We investigate in addition the effect of having an extra 
normalization factor /b for the 8 B neutrino flux, which is in this case assumed to be free. 
In what follows we will give explicit definitions of the \ 2 f° r the solar rates, the SK recoil- 
electron spectrum, the SK zenith angle distribution as well as the combination, to be used 
in our analysis. 



A. Calculation of the rates 

We calculate, for the two (three) generation framework, our theoretical predictions for 
the measured solar neutrino rates as a function of the two (three) mixing parameters for the 
gallium and chlorine experiments by folding the neutrino oscillation probability P2 g (-Pjg) 
with interaction cross section and the six solar neutrino fluxes corresponding to each reaction, 
pp, pep, 7 Be. 
Pinsonneault 
for simplicity. For gallium and chlorine experiments, we use the neutrino absorption cross 



8 B, 13 N and 15 as predicted by the standard solar model of Bahcall and 
U (BP98 SSM). Other minor sources such as 17 F and hep are not considered 



sections found in Ref. |K| and for the SK experiment, we use the new calculation of v e , 
scattering cross section on electrons which take into account radiative corrections j|6| . 



>,r 



In this case the expected event rate for Ga and CI detectors that should be compared to 
experimental data is : 



i4 heo = C x JdEa x (E)(P n (E)) IfsME) + / Be <M£) + E<M^) 



(16) 



where x=Ga,Cl, n=2g,3g, 4>b{E), <p Be (E) and <fij(E) with }=pp,pep, 13 ~N and 15 are the 
neutrino fluxes as a function of the neutrino energy taken from Ref. M, and C x is some 
normalization constant determined in such a way that _R^ heo is given in Solar Neutrino Unit 



(SNU). The numbers fe and fe e are respectively the normalization constants for the 8 B and 
the 7 Be neutrino fluxes, fe = fee = 1 corresponds to the BP98 SSM values. Here (P n ) means 
that the probability has been averaged on the neutrino path length to take into account the 
eccentricity of the Earth orbit around the Sun, after we have computed the probability as in 
Eq. (|7|), for two generations, or as in Eq. (|i3|) , for three generations, with the solar matter 
effect included. The integral above was performed starting at the energy threshold of each 
experiment. 

In our calculations, we have also included the effect due to the thermal broadening of 
about 1 keV of the 7 Be lines [47 when computing the capture rate by performing an extra 
average over the two 7 Be neutrino energy profiles given in Ref. |^7| . This is pertinent for 



large values of Amf 2 as discussed, for e.g. , in Refs. ||8|;|4T| . 

When fe is taken to be free we assume the power-law relationship between the 8 B and 
7 Be fluxes, given by Bc = (0 B ) 10 / 24 K7L i n order to renormalize the 7 Be flux by a factor 



fee as 

fee = /r 4 - (17) 

For the SK experiment, the expected solar neutrino event rate, normalized by the BP98 
SSM prediction, is given by, 

feJdE e JdE' e h(E' e ,E e )JdE„ME.) (^r( P n(E„)) +^i 1 ~ ( P n(^)) ]) 
Rt ^° = r ~ r ~ /._ , .„Aa v * ' 



J dE e J dE' e h{E' e , E e ) J dE u <f) B {E» 



dE' e 



(18) 

where E e is the observed recoil electron energy, E' e is the true recoil electron energy, h(E' e , E e ) 
is the electron energy resolution function taken from Ref. |4iJ and da u JdE' e , da v JdE' e are 
v e — e and z/ x — e (x = /i, r) scattering cross sections taken from Ref. |5(| . In the integral 
above we have used 5.5 MeV as the energy threshold for SK. 

The definition of the x? a tes function to be minimized is the same as the one used in 



Ref. [|50| which essentially follows the prescription given in Ref. |5T| . Explicitly, 

xLcs = E(^x heo - Rf>^{Rf co - Rf«l (19) 



xy 
x,y 



where x(y) runs through four solar neutrino experiments, Homestake, SAGE, 
GALLEX/GNO combined and Super-Kamiokande, the theoretical predictions R^m are 
given in Eqs. flTED and ( |TBD and the experimental values -R^y) are given in Table |. The 
error matrix cr xy contains the theoretical and experimental uncertainties according to Ref. 
BTj where theoretical uncertainties are taken from the ones given by the BP98 SSM ||. 



B. Calculation of the SK recoil-electron Spectrum 

For the spectrum shape analysis we first define the following quantity, 

/■£i max r r (da da 

/ dE e / dE' e h(E' eJ E e ) / dE„MEv) -JWF^ E »)) + ^ [ 1 " <?*{&>)) 



/" E i max r r da 

j Emin dE e J dE' e h(E' e} E e ) J dE v <h{E»)-^ 



(20) 



where _E. mm and _E. max are the minimum and maximum observed recoil electron energy in the 
i-th bin, starting at the 5.5 MeV threshold for the SK detector. In total we have 18 bins 
taken at 0.5 MeV intervals, except for the last bin which includes all the contribution above 
14 MeV. 

The definition of the Xspec function to be minimized is the same as the one used in 
Ref. |T7|, that is, 



X s 2 pec = E « he °-^ K K 2 K he0 -^ )' ( 21 ) 

ij=l,18 



where cry were computed as prescribed in Ref. |lj]] and Sf K are the experimental points 



whose numerical values are graphically reproduced from Ref. H. We will not use the new 
5 MeV point since the systematic errors of this point are still under study by the SK 
Collaboration. The extra normalization parameter a, which is always taken to be free in 
our analysis, is introduced because here we are only interested in fitting the shape of the 
spectrum. Moreover when we combine with the rates, it allows us to avoid double counting 
the information already taken into account in the rate analysis. 

For the SK spectrum analysis, as a good approximation, we simply assume that neutri- 
nos exit from the Sun as pure v e and use the vacuum probability formula in Eq. (^), for 
two generation, or in Eq. (0), for three generation, in this latter case removing from our 
expression the Sun matter effect. This is well justified by the discussion we presented in 
section III A, 



C. Calculation of the zenith angle dependence 

We define the Xzemth; f° r the zenith angle dependence, as follows, 

I p> vtheo T'obs \ 2 

2 _ v- iPfj -A ) (22) 

Azcnith ~~ Z^ 2 ' V / 

i=l,6 "Z,i 

where Z* heo is the i-th bin theoretical expectation, Z° hs is the i-th bin observed value, with 
5 night bins and 1 day bin which are graphically reproduced from Ref. |J and j3 is a free 
normalization constant introduced for the same reasons as in the case of the spectrum 
analysis. As the LVO solution does not imply in practice any zenith angle distortion of the 
data, this will simply provide an extra global constant that will increase the final combined 
Xmm- It is important to remark that this global constant will only affect the quality of the 
combined fit, bearing no influence on the computed allowed region. We obtained xLnith = 
4.3 with f3 = 0.47 for the zenith angle distribution. 

D. Combined Analysis 

Finally the combined x 2 to be minimized is simply defined as the sum of the individual 
ones: 

Acomb A-ratos ' A-spec ' Xzenith" V / 



E. Definitions of the confidence levels and A^of 

Using the y 2 functions defined in the previous subsections, for the two generation analysis 
where we have only two mixing parameters, (Am 2 2 , 6*12), we use the condition x 2 = x^ in + 
Ax 2 where Ax 2 = 4.61, 5.99 and 9.21 for 90, 95 and 99 % C.L., respectively, in order 
to determine the allowed parameter space. On the other hand, for the three generation 
analysis, in order to constrain the parameter space spanned by three variables (Am 2 2 , #12? 
#13), we determine the iso-confidence level surface by the condition x 2 — Xmin + ^X 2 where 
A X 2 = 6.25, 7.82 and 11.36 for 90, 95 and 99 % C.L., respectively, as in Ref. |35j. 

We note that in this work, we always determine the values of Xmin within the range 
10~ n eV 2 < Am 2 2 < 10 _8 eV 2 and therefore, do not take into consideration the region 



relevant for MSW solutions [[U|[T7[]. This means that in this work, "a priori 1 we assume 
that LVO is the solution to the solar neutrino problem, so adopting a different criteria than 
other authors, see for e.g. Ref. [36|, to draw the C.L. contours. 



Here, we also describe how we compute the number of degrees of freedom A^dof which 
is relevant to determine the goodness of fit (or C.L.). We compute A^dof as follows, 

-^DOF = ^"data — ^param — ^norm , (24) 

where Af data is the number of the data points we use in each analysis, N paxam is the number 
of mixing parameters to be constrained, two or three depending on the number of neutrino 
generations considered, and N norm is the number of extra free normalization factors we 
introduced in the analysis, i.e. fs in Eqs. ([16]) and (|18|), a in Eq. ( pT|) and (3 in Eq. (|22|). 



IV. RESULTS WITH THE TWO GENERATION SCHEME 

In this section we discuss our results for the analysis of the solar neutrino data in the 
context of two generations. 

A. Results with fixed /b 

We first present the case when we fix f B = 1, therefore f-Q e = 1. This corresponds to 
use the BP98 SSM flux values. To demonstrate the influence of the solar matter in the 
computed allowed regions, we will show here our results without and with the MSW effect 
in the Sun. 

In Figs. 0(a), (b) and (c) we show the allowed region in sin 2 # 12 — Am 2 2 parameter space, 
for the rates, spectrum and combined analysis, respectively, without taking into account 
any possible influence of the solar matter, hence all the plots are symmetric with respect to 
sin 2 #i2 = 0.5 because of the invariance of the pure vacuum probability under the transfor- 
mation 6*12 — >■ 7r/2 — #12. 

In Figs. |^(a) and (b) we show the allowed region in sin 2 Q 12 — Am\ 2 parameter space, 
for the rates and combined analysis, respectively, taking into account the MSW effect. As 
expected we see that in this case the shape of the allowed region is not symmetric with 
respect to sin 2 #12 = 0.5, for Am\ 2 > 3 x 10 -10 eV 2 , in good agreement with the results 
obtained in Ref. [II]. We also observe that there is a region allowed at 99 % C.L. for 



Am 2 2 > 4 x 10~ 10 eV 2 that only appears when matter effects are taken into account. 

No spectrum analysis with the solar matter effect was performed since we have checked 
that the results are essentially the same as the one presented in Fig. 0(b). This is because of 



the fact that nearly all neutrinos that are relevant to SK exit from the Sun as pure v e and 
so we can simply use the vacuum oscillation formulas as a very good approximation (see the 
discussions in section [II A|) . 



The values of the Xmin> the best fitted parameters, the number of degrees of freedom 
(-/V" D of), as well as the goodness of the fit, the C.L. in %, for the above discussed situations 
are shown in Table |TJ. In this table we also include this information for the local best fit 
point we found at Am 2 = 6.5 x 10~ n eV 2 and sin 2 # 12 = 0.2/0.8 which explain somewhat 
better the rates than the global best fit one. 

The MSW effect does not affect the values of the best fitted parameters for the spectrum 
and rates analysis, given two virtually symmetric solutions in sin 2 B\i. This symmetry is 
broken however when we combine rates and spectrum. Since the spectrum data prefers 
larger values of Am 2 2 , the combined best fit point will lay in an allowed island slightly 
deformed by the matter potential. 

Comparing the region allowed by the spectrum at 90 % C.L., given in Fig. || (b), with 
the one allowed by the rates at the same C.L., given in Fig. |3] (a), we can see that in general 
the region favored by the total rates are disfavored by the spectrum information. The SK 
spectrum data shows practically no energy dependent distortion, preferring either larger 
values of Am 2 2 , where the oscillating term gets averaged to one half, or smaller values of 
Am 2 2 , where the oscillating term vanishes. The total rates, however, prefer intermediate 
values of Am 2 2 , where the energy dependence of the solution provides the proper reduction 
to explain the four data points. Due to such disagreement between the total rates and the 
SK spectrum, we find it difficult to get a good fit for both at the same time. We note 
that the combined allowed regions are basically formed by the parameters which are more 
consistent with the spectrum data. 

In spite of this incompatibility we note that the quality of the combined Xmin/DOF = 
23.9/24 seems to be quite good. This can be understood in virtue of the fact that the 
oscillation parameters at this minimum can provide a very good explanation for the SK 
spectrum (x^n = 10.7) and zenith angle dependence (Xmin = 4.3), which are the majority of 
the statistical points, even though they give a generally poor explanation for the total rates 
(Xmm = 8.9). In fact the oscillation parameters at the combined minimum correctly predict 
the 71 Ga rate but for Homestake and SK the predicted values are substantially above what 
is measured by those experiments. 

B. Results with arbitrary /b 

For the case where we consider /b as a free parameter we also have performed the same 
analysis without and with matter effect taken into account. Although the allowed regions 
for the rates and the combined analysis increased a little bit they are substantially the same 
as the ones shown in Figs. |2| and [| so we do not show them here. 

Again, the values of the Xmin as we ^ as °f the best fitted parameters can be found in 
Table [II]. In comparison to the fixed /b case, the fit for the rates is substantially improved 
and the allowed parameter region became somewhat larger. The other qualitative features 
mentioned in the previous subsection remain unchanged. 

The combined fit also improved a bit by this extra freedom given x^/DOF = 21.7/23 
(see Table |ID, which seems to be even better than in the previous case. Allowing the 8 B 
flux to be free does not affect the fit of the SK spectrum and zenith dependence data, but 
it improves a little the fit for the rates (Xmm = 6-7) since we can now correctly predict the 



71 Ga and SK rates but still over estimate the Homestake value. 



V. RESULTS WITH THE THREE GENERATION SCHEME 

In this section we discuss our results for the analysis of the solar neutrino data in the 
context of three neutrino generations. Here we will always present the results taking into 
account the solar matter effect. 

A. Results with fixed /b 

As for the two generation case we first consider /b as fixed to be 1. In Fig. |] we show 
the regions, in sin 2 #12 — Am 2 2 parameter space, allowed by the rates for different values 
of sin 2 6*13. In each plot of Fig. [| we are presenting a cut of the allowed parameter space 
in three generations, that is, we show the allowed region in the sin 2 # 12 — Am 2 2 plane for a 
given value of sin 2 8 13 , which means that the Xmin va hie is common for all the plots shown. 
We found that for the rates the best fit occurs when sin 2 6 13 = 0. 

We note that, as in the case of the two generation analysis, the allowed regions that 
appear in the plots are asymmetric for Amf 2 ^3x 10~ 10 eV 2 because of the matter effect 
in the Sun. Moreover, there is a tendency that the matter influence will gradually start 
to be important at smaller values of Am 2 2 as sin 2 #13 increases. We can understand this 
by recalling that the effective potential here has been rescaled by cos 2 6*13 (see Eq. (|14D ), 
which means that as sin 2 #13 increases, the potential decreases and can start to be relatively 
relevant for lower values of Am 2 2 . 

In Fig. [5] the regions, in sin 2 9u — Am 2 2 parameter space, allowed by the spectrum for 
different values of sin 2 # 13 are displayed. Since we have used the vacuum oscillation formulas 
in our spectrum calculation, the allowed regions are symmetric with respect to sin 2 # 12 = 0.5. 
Again, as in the case of the rates analysis, we found the best fit at sin 2 #13 = 0. We note 
that although the best fit occurs when #13 = 0, the region excluded by the spectrum data 
become smaller as #13 increases. This can be qualitatively understood from the fact that 
the energy dependence of the probability becomes weaker as #13 increases (see Eq. (|13D) and 
consequently, the spectrum become flatter, which is consistent with the present SK data. 

In Fig. [^ we plot the combined allowed region for rates and spectrum, for fixed /b = 1. 
The values we have obtained for the best fitted parameters and Xmin as we ^ as the estimated 
C.L. of the fit are presented in Table |TJ. We have found that for the combined analysis, 
the best fit occurs when sin 2 6 , 13 = 0.12. Notwithstanding, a small difference in \ 2 such as 
Xmin(^i3 = 0) — Xmin(^i3 = 0.12) ~ 0.6 clearly bears no real statistical significance to the 
preferred non-zero #13. 

Finally we plot in Fig. |7| the spectrum calculated for the best fitted parameters for 
the rates, spectrum and combined analysis, for fixed values of sin 2 #13 , as well as the SK 
experimental data points. Here some comments are in order. One of the features of the total 
rate analysis is that several "local" best fit points, which are rather comparable in terms of 
X 2 values, exist. For 6 13 = we found the best fit at (Am 2 2 , sin 2 # 12 ) = (9.7xl0~ n eV 2 , 
0.35/0.65). We have noticed that as we increase the value of #13 these points do not always 
remain as best fit points but always remain as local best. Moreover, such two symmetric 
(with respect to 6*12 = 45°) best fit points move smoothly toward the direction of dyi = and 
#12 = 90°. The fitted spectrum curves for the rates, shown in Fig. [7], are, strictly speaking, 
indicating the "evolution" of the spectrum shape of these "best fit points" as 6^3 increases. 



We see that the best fit parameters for the rates do not produce a spectrum shape which 
is consistent with the SK data for lower values of sin 2 # 13 , but the combined fit, which is 
dominated by the spectrum weight in the x 2 is i n fairly good agreement with the data. 

As in the two generation case in spite of the incompatibility between rates and spectrum 
the quality of the combined fit where we found Xmin/DOF = 23.3/23 is quite good. Here 
again the oscillation parameters at the minimum can provide a very good explanation for 
the SK spectrum (x mill = H-l) and zenith angle dependence (x m in = 4.3), the majority 
of the statistical points, even though they give a generally poor explanation for the total 
rates (Xmin = 7.9), predicting total rates that are consistent with 71 Ga but inconsistent with 
Homestake and SK data at one or two sigma. 



B. Results with arbitrary /b 

The same plot as in Fig. || is presented for the combined analysis in Fig. |8|, when we 
considered /b to be free. We do not show the rates or the spectrum allowed regions in this 
case, since they are virtually the same as the ones shown in Figs. [| and |5|, respectively. 

We plot in Fig. £| \ 2 ~ X m in f° r the rates as a function of /b for three values of #i 3 . 
From this plot, we can see that as #13 becomes larger the allowed values of /b become more 
restricted, eventually, narrower than the range allowed by the BP98 SSM, 0.58 < /b < 1.57 
at 3 a [§]. The calculated spectrum for the best fitted parameters for the rates, spectrum 
and combined analysis, for various values of sin 2 # 13 , are similar to the case with /b = 1, 
except that the best fit point for the total rates can predict a less distorted spectrum shape 
at lower values of sin 2 #13 than in the previous case. We find it unnecessary to show this 
here. 

The combined fit improved a little with respect to the previous case given Xm m /DOF = 
21.7/22. Again allowing the 8 B flux to be free does not affect the fit of the SK spectrum 
and zenith dependence data, but it improves slightly the fit for the rates. 



C. Constraining #13 

The values of the Xmm as weu as °f the best fitted parameters for all the above cases are 
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shown in Table |TTTj. The combined x value includes Xzenith as f° r the two generation case. 
To illustrate the effect of the presence of the third neutrino in the \ 2 1 we plot i n Fig- P] 



the values of A% 2 = \ 2 ~ Xmin as a function of sin 2 #13 for the cases with /b fixed. From these 
plots, we see that there is a general tendency that the fit become worse as #13 increases, 
although we note that there is a local minimum for the fit with the rates and the combined 
data at sin 2 13 ~ 0.67. We also see that A% 2 for the rates increases more rapidly than that 
for the spectrum, which rises in fact quite smoothly. This can be understood as follows. 

As pointed out in Ref . |H| , naively, the fit for the rates is expected to become worse as #13 
increases since any energy dependence in the probability will become weaker (see Eq. flTB])) 
and thus, in general, the larger the value of #13, the similar will become the suppressions for 
all the solar neutrinos, leading to a gradually stronger inconsistency with the observed total 
rates. 

On the other hand, as we have already mentioned in the previous subsection |V A|, the 



loss of energy dependence will not compromise so much the fit for the spectrum data since 
the observed spectrum shape is consistent with a flat one. 



Finally, from Fig. 10, we can conclude that the solar neutrino data alone give the upper 
bound 9x3 < 63° at 95 % C.L. We have also performed the same analysis for arbitrary f B and 
obtained, in this case, a slightly weaker upper bound, 9i 3 < 65° at 95 % C.L. It is interesting 
to observe that this limit is quite similar to the one obtained in the case of a three flavor 



MSW solution to the SNP |35|| , although this MSW analysis was performed with different 
data. 



VI. PREDICTIONS FOR SEASONAL VARIATIONS 



In this section we discuss, in view of our previous presented results, possible seasonal 
variations that could be measured by the current Super-Kamiokande detector, as well as by 
the future Borexino |HJ and KamLAND [53j experiments, which will be sensitive to 7 Be 
neutrinos. It is clear that because of the eccentricity of the Earth orbit around the Sun, solar 
neutrino fluxes should vary as 1/L 2 as L varies with time. From Eq. (||) we can estimate 
that the flux variation in one year due to this effect is, 



A0, 



"h 



AeLc 



Ae<p u ~ 0.070, 



(25) 



which implies that the solar neutrino signals could vary as much as ~ 7 % in one year. 

We can see that if vacuum oscillation is assumed, "anomalous" time variation can be 
expected because the expressions of the probabilities in Eqs. (0),© and (0) also dependent 
on the distance L Hl4| , ^4|, |T5l . Such "anomalous" time variation is expected to be more 
prominent for 7 Be neutrinos |[55|j . This is because of the fact that 7 Be neutrinos are mono- 
energetic at E = 0.862 MeV (and also at E = 0.383 MeV, with much smaller flux) with 
AE ~ 1 keV [48], and the oscillation probability of such neutrinos can be very sensitive 
to the precise value of L, in contrast to the other neutrino fluxes such as pp and 8 B where 
the probability for these neutrinos must be averaged out over the neutrino energy. Detailed 
analyses on the seasonal variations for 7 Be neutrinos have been performed in Refs. p9|,|56|,|39 



In Fig. [jT| we plot the seasonal variation that is expected at SK for the best fitted 
parameters (global and local) of the two and three generation LVO solutions for E e > 5.5 
MeV jointly with the observed SK data f|. The global and local best fitted parameters 
correspond to the entries Combined and Combined (R), respectively, of Table 0, for the two 
generation, and of Table ITJ, for the three generation solution, with /b = 1 in both cases. 
SK presents the data without subtracting the expected "normal" seasonal variation due to 
the 1/L 2 dependence of the flux. For this reason we have included this extra effect in our 
curves in Fig. |TT[ We see that our predictions are currently consistent with the observed SK 
data. 

We also have computed the "anomalous" seasonal variation for the rate of 7 Be neutrino 
flux at 0.862 MeV that one can expect to be measured by the Borexino |5l| and the Kam- 
LAND [53] experiments using the best fitted parameters (global and local) for the LVO 
solution with two and three flavors. Here we have subtracted the expected "normal" sea- 
sonal variation due to the 1/L 2 dependence of the flux. We note that the predictions for 
Borexino and KamLAND we have obtained are visually indistinguishable, and therefore, 



we only plot the case of Borexino in Fig. [12]. In this figure our predictions are normalized 
to the BP98 SSM 7 Be neutrino flux value, i.e. /b c = 1- If /Be is different from unity, our 
predictions must be renormalized accordingly, keeping the same shape. As we can see from 
these plots, the best fit for two or three neutrino generations give similar shapes, also large 



time variation is obtained, which should be observable at both of these experiments. Such 
drastic variation is due to the fact that the oscillation wavelength for 7 Be neutrinos for our 
global best fitted values of Am 2 12 = 4.6 x 10" 10 eV 2 and 6.6 x KT 10 eV 2 , is about 3-4 % of 
the mean Sun- Earth distance (see Eq. (|4])), which is comparable to the Sun-Earth distance 
variation due to the orbit eccentricity. We have taken into account in our calculations the 
finite width of the 7 Be line IE7 . 



Such "anomalous" seasonal variation can be a clear signature of the vacuum oscillation, 
which does not depend on any detail of the SSM nor on unknown experimental systematic 
errors 



U 



VII. ANALYSIS WITHOUT CHLORINE DATA 



Finally, in this section, following Refs. 5e , we have further investigated the impact of 



removing from our analysis the chlorine data, since Homestake is the only radiochemical 
experiment which has not been calibrated with a radioactive source. This is certainly an 
extreme case but could be worthwhile to be discussed. 



A. Two generation case 

In Figs. [H|(a) and (b) we show our results for the rates alone and for the combined 
analysis, respectively, in the case of two generations. We observe here that the rates alone 
allow for a broader range of solutions in the sin 2 #12 — Am 2 2 plane, compared to Fig. |^(a) 
where the chlorine data was included, making it more consistent with the SK spectrum 
measurement. The combined analysis without chlorine data reflects this fact by allowing for 
larger values of Am 2 2 even at 90 % CL. 

The values of the Xmin as wen as °f the best fitted parameters here are shown in Tables 
rV] and Table [V]. Again the combined x 2 value also includes xLnith- We remark that without 



chlorine data the quality of the combined fit improves since there are significant overlap 
between the allowed parameter region which give good fit to the total rates and SK spectrum. 
We can see this by comparing Figs. |2](b) and |l3|(a). This is in contrast to the case with 
chlorine data. 



B. Three generation case 

We repeat the same procedure for the three generation case. In Figs, in] we show our 



results for the combined three flavor oscillation analysis, for several values of sin 2 # 13 . Again 
we see that if one ignores Homestake's data completely a larger range of values of Am 2 2 and 
sin 2 #12 can provide a good explanation for the measured rates as well as for all the solar 
data combined. 

In Fig. [IUJ we also show the values of A% 2 = x 2 ~ Xmin as a function of sin 2 # 13 for the 
analysis without 37 C1 data. This curve has the same features of the one for the case where 
all the data was considered in the study, and allows us to put an upper limit on #13 to be 
less than ~ 66° at 95 % C.L.. 



VIII. DISCUSSIONS AND CONCLUSIONS 

We have re-examined the status of the SNP in the light of the LVO solution, using 
the most recent solar neutrino data as well as the predictions of the BP98 SSM. We have 
analyzed the solar neutrino data in the context of two and three neutrino flavor oscillations 
in vacuum, extending previous analyses to covered the range 10~ n < Am 2 2 < 10~ 8 eV 2 . 
In doing this we have included the MSW effect in the Sun, which plays a relevant role in 
the calculations when Am 2 2 /E > few xlO -10 eV 2 /MeV. When we have considered three 
neutrino generations we assumed a mass hierarchy such that Am 2 3 should be large enough 
to be relevant for the atmospheric neutrino problem. 

We have found that the MSW effect in the Sun significantly modify the allowed parameter 
space for Am 2 2 larger than ~ 3 x 10~ 10 eV 2 in agreement with the result obtained in Ref. JIT . 
We, however, found that values of the best fitted parameters are not practically affected by 
the presence of the solar matter effect (see Table II). 

In the two generation LVO solution we found that the rates prefer lower values of Am 2 , 
i.e. Am 2 2 < 3 x 1CT 10 eV 2 at 90 % C.L., while the SK electron recoil spectrum data excludes 
such favored parameter region and prefers Am\ 2 > 6 x 10~ 10 eV 2 . In the combined analysis 
the weight of the spectrum prevails. 

The disagreement between the best fitted rates and the spectrum predictions for the LVO 
solution in the case of two generation is, nevertheless, less striking when three neutrinos are 
considered. When #i 3 is small, the picture described above applies but as 0i 3 increase the 
region favored by the total rates are less excluded by the spectrum data compared to the 
two generation case (see for e.g. , Fig. |] and Fig. |5| for the case sin 2 # 13 = 0.15). When /b 
is arbitrary there is even better agreement between spectrum and rates, the regions allowed 
by them already superpose for sin 2 #13 = 0. 

Judging by the values of the xLm we have obtained in our combined analyses one could 
conclude that the LVO solution is in quite good shape. We point out, however, that this fact 
should be understood as follows: LVO can provide a very good explanation for the spectrum 
and zenith angle dependence in SK as well as for the total rate of the 71 Ga experiments, 
these are the majority of the statistical points used in our analyses, but provides a poor 
explanation for the total rates at Homestake and SK. The SK rate can also be well explained 
if one allows for /b ~ 0.74. 

We found that the fit to the total rates as well as SK spectrum are not essentially 
improved due to the possibility of oscillating into a third neutrino (see Fig. |j~0[) . In particular, 
we obtained the upper bound on 6*13 to be less than ~ 63° at 95 % C.L. by the solar neutrino 
data alone. We should remark that this bound is substantially weaker than the bound 
obtained by combining the Super-Kamiokande atmospheric neutrino observations with the 
CHOOZ reactor experiment limit |23,25|. Hopefully future neutrino oscillation experiments 
at neutrino factories will be able to give a much more precise information on # 13 



Note added: While we were completing this work, we became aware of two similar works 
that appeared in hep-ph [|58| , |5971 . 
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TABLES 

TABLE I. Solar neutrino rates observed by Homestake, SAGE and GALLEX/GNO combined 
as well as theoretical predictions for the the Standard Solar Model by Bahcall and Pinsonneault 
P] . For SK we show the ratio of the observed flux over the prediction of the BP98 SSM. 



Experiment 



Observed Rate Ref. BP98 SSM Predictions f| 



Units 



-1.0 

129Jl| 



Homestake 

SAGE 

GALLEX/GNO 

Super-Kamiokande 



2.56 ± 0.23 

75.4 ± 7.6 

74.1 ± 6.8 

0.465 ± 0.015 



7.7; 



129 



4-8 



1.00 



+0.19 
-0.14 



SNU 
SNU 
SNU 
5.15 x 10 6 cm- 2 s 



TABLE II. The best fitted parameters and x m in as wen as C.L. (in %) for the 2 generation 
LVO solution to the SNP. The number of degrees of freedom (./Vdof) are also indicated. The 
matter effect was taken into account unless it is indicated in parentheses. Note that for spectrum 



analysis matter effect was neglected, as it is a good approximation (see section [II A] , [III B|) . We 
also present the values of the local best fit which explain better the total rates in the line indicated 
as Combined (R). 



Case 


Am 2 2 x 10 10 eV 2 


sin 2 6»i2 


/b 


Y 2 - 

A, mm 


-^DOF 


C.L. (%) 


Rates 


0.97 


0.35/0.65 


1 (fixed) 


0.25 


2 


88.2 


Spectrum 


6.5 


0.38/0.62 


- 


10.5 


15 


78.7 


Combined (w/o matter) 


6.6 


0.35/0.65 


1 (fixed) 


25.1 


24 


40.0 


Combined 


6.6 


0.64 


1 (fixed) 


23.9 


24 


46.7 


Combined (R) 


0.65 


0.2/0.8 


1 (fixed) 


28.0 


24 


26.0 


Rates 


0.88 


0.34/0.66 


1.21 


0.01 


1 


91.6 


Combined (w/o matter) 


6.6 


0.31/0.69 


0.72 


22.7 


23 


47.8 


Combined 


6.6 


0.66 


0.74 


21.7 


23 


53.8 


Combined (R) 


0.65 


0.2/0.8 


0.86 


27.5 


23 


23.5 



TABLE III. The best fitted parameters and Xmm f° r the 3 generation LVO solution to the SNP. 



Case 



sin 013 



0i3 Am? 9 x 10 10 eV 2 



'12 



sin 



7 12 



y 2 - 

A-min 



iVi 



DOF 



C.L. (%) 



Rates 


0.0 


0.97 


0.35/0.65 


1 (fixed) 


0.25 


1 


61.7 


Spectrum 


0.0 


6.5 


0.38/0.62 


- 


10.5 


14 


72.5 


Combined 


0.12 


4.6 


0.75 


1 (fixed) 


23.3 


23 


44.3 


Combined (R) 


0.14 


0.66 


0.15/0.85 


1 (fixed) 


26.1 


23 


29.6 



Rates 


0.0 


0.88 


0.34/0.66 


1.21 


0.01 





- 


Combined 


0.0 


6.6 


0.66 


0.74 


21.7 


22 


47.7 


Combined (R) 


0.1 


0.65 


0.14/0.86 


0.82 


25.5 


22 


27.4 



TABLE IV. The best fitted parameters and x m in f° r ^ ne 2 generation LVO solution to the SNP 
without the chlorine data. 


Case 


Am 2 2 xl0 10 eV 2 sin 2 12 f B X mi n 


N DOF C.L. (%) 


Rates 
Combined 


0.97 0.33/0.67 1 (fixed) 0.01 
6.5 0.69 1 (fixed) 18.0 


1 
23 


92.0 

75.7 


TABLE V. The best fitted parameters and Xmin f° r the 3 generation LVO solution to the SNP 
without the chlorine data. 


Case sin 2 13 Am 2 2 x 10 10 eV 2 sin 2 12 f B 


y 2 - 

Amin 


iV D OF C.L. (%) 


Rates 
Combined 


0.0 0.97 0.33/0.67 1 (fixed) 
0.1 4.6 0.76 1 (fixed) 


0.01 

17.4 



22 74.1 



FIGURES 

FIG. 1. Contours of iso-survival probability of v e at the surface of the Sun as a function of 
Amf 2 /E and sin 2 B\2 (a) considering pure vacuum oscillation, and (b) when matter is taken into 
account. For the case (b) we numerically integrated the evolution Eq. (|I|) assuming that neutrinos 
are created at the center of the Sun. 

FIG. 2. Region of (sin 2 12 , Am 2 2 ) allowed by (a) the total rates from 37 C1, 71 Ga and SK 
experiments, (b) the SK spectrum and (c) the combined analysis of rates + SK spectrum in the 
long-wavelength vacuum oscillation scenario for 2 neutrino flavors. In the case of the SK spectrum 
(b) it is the inner part of the contours which is excluded by the data. The best fit points are shown 
as black circles. Solar matter effect is not taken into account. 



FIG. 3. Same as in Figs. || (a) and (c) but with the solar matter effect. 

FIG. 4. Region of (sin 2 #12, Am 2 2 ) allowed by the total rates from 37 C1, Ga and SK for 
various values of sin 2 #13 for the long-wavelength vacuum oscillation solution to the SNP with 3 
neutrino flavors. 

FIG. 5. Region of (sin 2 #i 2 , Am 2 2 ) allowed by the SK spectrum for various values of sin 2 #13 
for the LVO solution to the SNP with 3 neutrino flavors. It is the inner part of the contours which 
is excluded by the data. 

FIG. 6. Same as in Fig. || but for the combined analysis of total rates + SK spectrum. 

FIG. 7. Expected SK spectrum using the best fitted values of (sin 2 9\2, Am 2 2 ). We note that 
the spectrum curves determined only by the rates are adjusted in such a way that x 2 defined in 
Eq. ( pi|) takes minimum values after we determine the values of (sin 2 612, Am 2 2 ) by the fit only 
with the rates. 

FIG. 8. Same as in Fig. y but with arbitrary 8 B neutrino flux normalization /b- 

FIG. 9. x 2 ~ Xmin i s plotted as a function of /b for various values of sin 2 #13 . We also show, 
by the vertical lines, the range of allowed values of /b at one and three standard deviations as 
predicted by the BP98 SSM |8|. 



FIG. 10. x 2 ~ Xmin as a function of sin 2 #13, for the case when /b is fixed to be equal to 1 for 
the rates (R), spectrum (S), combined analysis (R+S) and for the combined analysis without 37 C1 
data (R+S noCl). 

FIG. 11. Expected seasonal variation at SK for the best fitted parameters of the two and three 
generation LVO solutions to the SNP. We also show in the plot the experimental data points. We 
have not subtracted the effect of the "normal" seasonal variation (due to ~ 1/L 2 dependence) 
expected in the absence of any oscillation. 

FIG. 12. Same as in Fig. 11 but for Borexino and KamLAND. Here we have subtracted the 
"normal" seasonal variation due to the 1/L 2 dependence of the flux. 



FIG. 13. Region of (sin 2 #12, Ara^) allowed by (a) the total rates from 71 Ga and SK, (b) the 
combined analysis of rates + SK spectrum in the LVO scenario for 2 neutrino flavors. The best fit 
points are shown as a dark circle. Here we have ignored the information from CI. 

FIG. 14. Region of (sin #12, Am 2 2 ) allowed by the combined analysis of rates + SK spectrum 
for various values of sin 2 #13 for the LVO solution to the SNP with 3 neutrino flavors. Here we 
have ignored the information from 37 C1. 
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